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Abstract 


The  ablation  characteristics  of  a  thick  walled  spherical  shell 
considering  radial  heat  flow  is  investigated.  The  outer  sur¬ 
face  of  the  shell  is  subjected  to  a  time-dependent,  point- 
symmetric,  radial  heat  input  and  the  inner  surface  is  insulated. 
The  melt6d  material  is  assumed  to  be  immediately  removed  upon 
formation. 

Two  approximate  solutions  to  the  problem  are  developed  and  com¬ 
pared.  The  first  method  is  an  adoption  of  Citron's  solution 
[l]  which  is  based  on  the  assumption  that  the  temperature  dis¬ 
tribution  through  the  thickness  of  the  body  may  be  expressed 
in  a  Taylor  series  expansion  about  the  melting  surface  at  any 
time.  The  second  method  employs  the  heat  balance  technique 
suggested  by  Goodman  [2]  which  satisfies  the  heat  conduction 
equation  on  the  average. 

Numerical  calculations  are  performed  for  an  aluminum  sphere 
having  a  7-3/4  inch  outside  diameter  and  a  one-inch  wall  thick¬ 
ness,  subjected  to  a  constant  heat  flux.  Constant  thermal 
properties  are  assumed  for  these  calculations. 


Numbers  inside  the  square  brackets  refer  to  the  References. 
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Symbols 


a 

b 

B* 

c(T) 

c 

G 

h 

k(T) 

k 

L 

* 

M 

Q(t) 

Q 

CL 


rs(t) 


initial  outside  radius  of  spherical  shell 


t 

Si 

T(r,t) 

T0(r,t) 


inside  radius  of  spherical  shell 

Q„[s  -  b] 


non-dimensional  parameter; 


specific  heat  of  material,  a  function  of  temperature 

c(T) 

non-dimensional  specific  heat  parameter;  77-7 — r 

rs 

temperature  integral;  j  u  dr  , 

b 

surface  heat  transfer  coefficient 

thermal  conductivity  of  material,  a  function  of 
temperature 


non-dimensional  thermal  conductivity  parameter; 
latent  heat  of  fusion 


k(T) 

k(Tm> 


non-dimensional  parameter; 


e(TB)[T  -Tt] 


Q(t) 


heat  flux  on  the  outer  surface  of  body  [for  aerodynamic 
heating  Q( t) *  h(T#  -  Tw)] 

non-dimensional  heat  flux  parameter; 

heat  flux  at  t  =  t 

m 

radial  distance  # 

outside  radius  of  the  spherical  shell  at  any  time 

a  -  rs 

non-dimensional  outside  radius  parameter;  ~  7 “ * 
time 

melt  time,  time  at  which  the  melting  temperature  is 
first  reached  on  the  outside  surface  ( r »  a ) 

temperature  (t > t  ) 

pre%-melt  temperature  (t£t  ) 

ii 


If 

p 

T 


constant  initial  temperature  of  the  body 

melting  temperature  of  the  material 

stagnation  temperature  of  the  external  flow  field 

■ 

outside  surface  temperature  [for  t£.tm,  S  Tffl] 

temperature  transformation;  u=tr(uo-TQr) 

non-dimensional  space  transformation;  r  -  b/rg  -  b 

non-dimensional  t’emperature  parameter;  T-T,/Tm-T^ 

non-dimensional  surface  heat  transfer  coefficient 

[*o &  haA] 

thermal  diffusivity;  k/pc 


density 


non-dimensional  time  parameter; 


iii 


Introduction 


Upon  re-entry  the  extremely  high  heat  inputs  have  major  effects 
on  the  structural  integrity  of  the  space  vehicle.  The  heat  in¬ 
puts  are  of  such  magnitude  that  the  melting  temperature  of  the 
space  vehicle  may  be  reached  on  the  outside  surface.  Melting  of 
the  surface  will  then  occur,  and  the  aerodynamic  shearing  forces 
will  tend  to  remove  this  molten  material  from  the  original  struc¬ 
ture.  This  destructive  process,  which  is  generally  termed  ab¬ 
lation,  has  been  used  quite  successfully  for  beneficial  purposes, 
such  as  the  heat  shield  design.  This  technique  employs  coating 
the  space  vehicle  with  a  shielding  material.  During  re-entry 
the  ablation  of  this  shielding  material  absorbs  a  great  per¬ 
centage  of  the  heat  input  to  the  entire  body  and  thereby  protects 
the  load-carrying  structure  and  interior  of  the  vehicle  from  heat 
damage.  It  is  of  importance  to  the  designer  to  obtain  an  accu¬ 
rate  prediction  of  the  amount  of  shielding  material  needed  to  ab¬ 
sorb  a  given  heat  input  and  to  determine  the  corresponding  tem¬ 
peratures  and  thermal  stresses  within  the  shield  and  back-up 
structure . 

The  treatment  of  any  problem  in  heat  conduction  involving  more 
than  one  space  variable  is  a  complicated  procedure.  With  the  in¬ 
clusion  of  ablation  of  the  outside  surface,  the  problem  becomes 
non-linear,  even  in  the  simplest  cases.  For  this  reason,  all  the 
available  literature  investigated  treated  problems  in  one  space 
variable. 

The  basic  modern  work  on  ablation  was  carried  out  in  1948  by 
Landau  [3]  in  which  he  set  up  the  equations  and  boundary  condi¬ 
tions  for  ablation  of  a  finite  and  semi-infinite  slab  and  carried 
out  the  solution  for  the  semi-infinite  slab  under  constant  heat 
input  and  thermal  properties  with  the  use  of  computing  machinery. 
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It  was  assumed  that  the  melted  material  was  immediately  removed 
upon  formation. 

Lotkin  [4]  developed  a  numerical  procedure  for  the  solution  of 
an  ablating  finite  slab.  Immediate  removal  of  the  melt  was  as¬ 
sumed,  as  were  a  time-dependent  heat  flux,  and  temperature- 
dependent  thermal  properties.  Again  high  speed  computing  machines 
were  essential  to  the  solution  of  the  problem. 

A  computer  technique  for  the  solution  of  the  slab  problem  under  a 
time-dependent  heat  input  was  developed  by  Ehrlich  £5].  In  this 
case,  however,  the  melted  material  was  assumed  to  remain  in  con¬ 
tact  with  the  original  body,  necessitating  satisfying  the  heat 
conduction  equation  in  both  the  melted  and  unmelted  material. 

Dewey,  Schlesinger  and  Sashkin  [6]  developed  a  numerical  solution 
for  a  cylinder  of  finite  thickness  under  radial  heat  input  with 
one  moving  boundary  and  variable  thermal  conductivity.  A  simple 
mathematical  model  for  aerodynamic  ablation  was  developed  by 
Goodman  [7],  and  Adams  [8]  discussed  the  important  ablation  para¬ 
meters  and  the  various  solutions  which  have  been  obtained,  and 
compared  these  to  experimental  results.  Excellent  background 
material  for  the  entire  heat  conduction  problem  as  well  as  some 
specific  solutions  related  to  changes  of  phase  can  be  found  in 
Car.slaw  and  Jaeger  [9]. 

In  an  attempt  to  arrive  at  solutions  of  the  ablation  problem 
without  the  use  of  high-speed  computing  machinery,  a  few  approxi¬ 
mate  techniques  have  been  developed.  Goodman  [2]  employed  the 
heat  balance  technique  for  the  solution,  in  closed,  analytical 
form,  of  the  ablating  s em i -in f in i te  slab  and  compared  these  re¬ 
sults  with  Landau's  exact  solution.  Citron  [1,10]  developed  two 
techniques  for  the  solution  of  the  finite  slab.  One  method  [1] 
was  based  on  the  assumption  that  the  temperature  distribution 
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through  the  thickness  of  the  slab  could  be  expressed  in  a  Taylor 
series  expansion  in  space  about  the  melting  surface.  This 
yielded  an  ordinary,  non-linear,  differential  equation  in  terms 
of  the  melt  depth  as  a  function  of  time.  The  solution,  which 
must  be  obtained  numerically,  could  be  readily  computed  on  a  desk 
calculator.  Citron's  second  method  [10]  consisted  of  reducing 
the  n<jn-linear,  partial  differential  equation  into  two  ordinary 
differential  equations,  one  linear  and  one  non-linear,  in  terms 
of  a  temperature  function  and  the  melt  depth.  The  solution  could 
be  obtained  by  successive  approximations  utilizing  a  desk  calcu¬ 
lator.  Boley  [ll]  cons  idered  •  st i  1 1  another  procedure  for  the 
slab  with  immediate  removal  of  the  melt  and  constant  thermal 
properties.  An  ordinary  int egro-d i f f er ent ia 1  problem  which  can 
be  solved  numerically  'or  in  series  form  for  the  exact  solution  of 
the  melting  problem  was  developed.  A  solution  for  the  semi¬ 
infinite  slab  under  constant  heat  input  was  obtained. 

Interaction  of  the  ablating  material  with  the  external  flow  field 
was  also  considered  by  certain  investigators.  When  a  material 
ablates,  the  molten  material,  or  the  vaporized  material  in  the 
case  of  a  subliming  solid,  or  a  mixture  of  vapor  and  liquid,  is 
injected  from  the  body  into'the  surroundings  of  the  body.  In  the 
case  of  re-entry,  the  injection  is  made  into  the  boundary  layer 
and  this  produces  the  added  beneficial  effect  of  reducing  the 
heat  input  to  the  body.  Swann  and  South  [12],- Lew  and  Fanucci 
[13],  Fleddermann  and  Hurwicz  [14],  and  Sutton  [15],  include  in¬ 
vestigations  of  such  phenomena.  In  Ref.  [16],  Economos  includes 
this  factor  in  calculating  ablation  of  semi-infinite  slabs  of 
plastic  materials,  6dch  as  Lucite,  and  compares  the  results  with 
experimental  data. 

In  Ihis  report,  an  attempt  has  been  made  to  determine  the  solution 
of  the  ablation  problem  for  a  thick  walled  hollow  sphere,  using 


the  Goodman  [2]  and  Citron  [1]  techniques. 

The  two  methods  have  been  compared  in  an  attempt  to  evaluate  the 
inherent  attributes  of  each  and  to  determine,  if  possible,  which 
method  is  more  practical  in  actual  application.  It  was  assumed 
that  the  outer  surface  of  the  sphere  was  subjected  to  a  pre¬ 
scribed  point-symmetric,  time-dependent,  radial  heat  flux  and 
that  the  inner  surface  of  the  sphere  was  insulated.  The  initial 
temperature  distribution  in  the  body  was  constant  and  lower  than 
the  melting  temperature  of  the  material.  Thermal  properties 
were  considered  to  be  functions  of  temperature  (for  the  Citron 
method)  and  it  was  also  assumed  that  the  melted  material  was  im¬ 
mediately  removed  upon  formation. 

Numerical  calculations  of  melt  depth  and  temperature  distribu¬ 
tions  as  functions  of  time  were  carried  out  for  an  aluminum 
sphere  of  1"  wall  thickness  and  7-3/4"  outside  diameter  under 
aerodynamic  heat  inputs  corresponding  to  hypersonic  flow.  Materi¬ 
al  properties  were  assumed  constant  for  these  calculations. 
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Theoretical  Analysis 

The  problem  treated  is  the  ablation  of  a  hollow  sphere  initially 
at  a  uniform  temperature  (T^)  which  is  lower  than  the  melt  tem¬ 
perature.  It  is  assumed  that  the  outer  surface  (r=a)  is  subjec¬ 
ted  to  a  time-dependent,  point-symmetric,  radial  heat  input  Q(t), 
that  the  inner  surface  (r=b)  is  insulated,  and  that  the  material 
has  temperature-dependent  thermal  properties.  The  molten  ma¬ 
terial  is  assumed  to  be  immediately  removed  so  that  the  outer 
boundary  of  the  body  is  always  considered  to  be  at  the  melt  tem¬ 
perature. 

The  overall  analysis  can  be  divided  into  two  separate  investiga¬ 
tions!  i)  the  pre-melt  analysis,  and  ii)  the  melt  analysis* 


i)  Pre-melt  Analysis! 

This  analysis  covers  the  time  interval  between  the  initial  time 
(t*0),  when  the  body  is  at  a  constant  initial  temperature  (T^), 
and  the  melt  time  (t-tffl),  the  time  at  which  the  melting  tempera¬ 
ture  (T  )  of  the  material  is  first  reached  at  the  heated  surface, 
m 

Attainment  of  this  solution  is  necessary  for  the  melt  analysis, 
for  in  order  to  continue  the  solution  of  the  problem  once  melting 
begins  it  is  necessary  to  know  the  time  ( t  )  when  melting  starts 
and  the  temperature  distribution  through  the  thickness  of  the 
body  at  this  time. 

Before  melt  temperature  is  reached  on  the  outer  surface  of  the 
hollow  sphere,  the  heat  conduction  equation  isi 

8To  19  2  ®Io 

Pc  ~  =  72^7(kr  1  b<r<*  U) 
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under  the  conditions-. 


a  ) 


b) 


c ) 


Vr’0)=Ti 

3T,, 

Q(t)  =  k(-rr) 

8r  a  »  t 


(2) 


Many  solutions  of  the  pre-melt  problem  exist  (see,  e.g.,  Refs. 
[9]  and  [17]),  and  thus  in  general  the  melt  time  and  the  corres¬ 
ponding  temperature  distribution  can  be  readily  obtained. 


ii)  Melt  Analysis; 

When  melting  begins  it  is  assumed  that  the  outer  surface  of  the 
sphere  always  remains  at  the  melt  temperature  and  that  the  melt 
is  immediately  removed.  The  total  heat  input  to  the  body  is  now 
divided  into  two  parts.  One  portion  enters  the  solid  while  the 
other  accounts  for  the  latent  heat  of  fusion  absorbed  in  the  ab¬ 
lation  process. 

The  heat  conduction  equation  is: 

9T  1  a  ,  2^! 

pc  at  =  7?  a7Ur  a7>  5  b<r<rs(t)  U) 

where  rg(t)  is  the  varying  outside  ra'dius  of  the  sohere. 


The  conditions  are: 


-  7  - 


a) 

T(r*tm) 

=  To(r’t, 

b) 

T(rs,t) 

s  Tm 

c) 

Q(t)  = 

,3T 

d) 

8r  b,t 

=  0 

e). 

rs(tm) 

=  a 

(4) 


An  auxiliary  condition  can  be  imposed  from  a  consideration  of 
continuity  of  heat  input  at  time  t=  t^.  From  the  pre-melt  ana¬ 
lysis  the  boundary  condition  on  the  heat  input  at  t *  t  from 
Eq.  (2b)  isi 


8T, 

Q(t  )  =  k(T  )(-T— ) 
m  m  o  r 


a  »tn 


and  from  the  melt  analysis  the  Boundary  Condition  from  Eq.  (4c) 
is  i 


Now  , 
ity  t 


AT  dr. 

„  ,  -  PL<-ar’  ' 

VV’V  *m 

since  r  (t  )  s  a,  and  T(  r»  t  )  *  T  (  r ,  t_  ) »  then  for  continu- 
s  m  mom 


The  solution 
(  5  )*  wi  1 1  now 
techniques . 


drs 

m 


(V 


of  the  problem,  described  by  Eqs.  (3),  (4),  a.nd 
be  treated  by  two  different  approximate  numerical 
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Method  It 

In  this  approach,  the  technique  utilized  by  Citron  [l]  for  a  slab 
is  applied  to  the  spherical  shell.  Material  properties  are  as¬ 
sumed  to  be  functions  of  temperature.  The  method  consists  of  ap¬ 
plying  a  transformation  which  allows  the  consideration  of  a  body 
of  constant  unit  thickness  at  all  times  in  lieu  of  a  body  of 
varying  thickness,  and  then  e\pressing  the  temperature  distribu¬ 
tion  at  any  time  in  this  unit  body  by  a  Taylor  series  expansion 
in  space  about  the  melting  surface. 

For  the  spherical  shell  the  transformation  used  is 


so  that  at 
and  at 


Z  = 


r  -  b 
r  -  b 


r=  V 


2-  1 
r =  b  ,  Z=  0 


(6) 


Using  this  transformation  along  with  the  following  non- 
dimensional  parameters 


«(Tm)(t  ”  tm ) 

T  =  - 5 - 

(a  -b)2 


Rs(t) 


T  -Tt 

e( z,t)  *  ^ — ry- 

m  i 


(?) 


k 


k(T) 


-  c(T) 

c  =  C<V 


where  Qq  is  the  heat  input  at  tm»  The  heat  conduction  Equation 
(3)  becomes 


i# 
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c,  ,2  9©  c  •  8$ 

£  ( 1  “Rs)  at  +  1  "Rs}  Rs  8Z 


2(1  -Rs) 


o-",’2  +  wi 


96 

az 


I  ii  /®i\' 

^  de  (az' 


(8) 


where  (  )  denotes  differentiation  with  respect  to  T,  and  the 
conditions  (4a  to  4e)  become 


a ) 

8(Z,t)  = 

6q(Z,t)  at 

T  = 

b) 

6(Z,t)  = 

1  at 

z  = 

c) 

(fii)  = 

92  1,T 

(l  -RJC-T  * 

S  B 

• 

Rs 

— ] 

*J 

M 

d) 

(2£)  = 
QZ  0,T 

0 

e) 

Rs(0)  = 

0 

drs 

while  Eq.  (5),  (j£~)  =  0,  becomes  Rg(0)  =  0 

*m 


and 


.  MT„)[Tm-Ti] 
*  0oCa-b] 

Tm )  CTm  “  T  ^  ] 
m  -  : 


(9) 


It  is  now  assumed  that  0(Z»t)  can  be  expressed  as  a  Taylor 
series  expansion  in  space  about  the  melting  face  Z=  1.  That  is, 


e(  z,t)  =  e(  i  ,-r)  +  (|-) 


1  *  T 


(Z-l)  +  (^-|) 

az  i,t 


(z  - 1)2 


21 


(10) 
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The  first  term  and  the  coefficient  of  the  second  term  are  known 
from  conditions  (9b'and  9c) 


6(1, t)  =  1 


(fii) 


1  »T 


=  (.1  -*R.)[-T 


The  coefficient  of  the  third  term  can  readily  be  found  by  evalu¬ 
ating  Eq.  (8)  at  Z  =  1.  The*result  is 


(^-|)  =  (*1  -R8)2f[-T  “  ~ 

3Z *  1 1 T  *  B  M  8 


1  —  R  4-  !?/ 9 

1  Rs  +  l.-b/a 


dk 

de 


) 

l.T 


} 


By  successively  differentiating  .Eq.  (8)  with  respect  to  Z, 

a3e  a4e  ane 

•a  »  a  »  •••' i  . ""  »  can  be  obtained  and  these  can  be  evalu- 

azJ  az  azn 

ated  at  (1,t)  to  form  the  remaining  n-2  coefficients  of  thfe 
Taylor  series  expansion.  The  result  is  that  an  expression  for 
0(Z,t)  is  obta ined 'conta ining  Z  and  and  its  first  K  deriva¬ 
tives  when  2K  or  2K  +  1  terms  of  the  Taylor  series  are  included. 


Now  if  condition  (9d),  which  is 


(az)  =  0 


0,T 


is  applied  to  Eq.  ( 1£ ) ,  then 


„  .as.,  ,aV  •  i  ,a3e 

°  ‘  (82>  '  ♦  2T<-3>, 

i»t  az  i*t  az  i,t 


( 11 ) 


11 


-3$ 


and  since  (4—)  }  ( 

82  l,x 


82^  1,t 


8ne, 


;  (“) 


9Z  1  *t 


have  all  been  ex¬ 


pressed  in  terms  of  R  .  R  ,  etc.,  the  final  result  is  that  a 

s  s 

non-linear,  ordinary  differential  equation  involving  R?  and  its 
derivatives  alone  has  been  obtained.  It  should  be  noted  that 
this  equation  will  always  be  linear  in  the  highest  order  deri¬ 
vative  term  (when  more  than  three  terms  are  used  in  the  expan¬ 
sion)  . 


The  series  expansion  of  the  temperature  is  terminated  after  2K 

or  2K  +  1  terms,  thus  leading  to  a  differential  equation  of  Kth 

order.  Thus  in  addition  to  the  initial  conditions  R  (0) =  0, 

•  ® 
and  R  (0)  =  0,  it  is  required  to  obtain  K-2  additional  condi- 

tions.  These  K-2  values,  R  (0),  R  (0)  ....,  R  (0)  are  obtained 

s  s  s 

by  matching  the  initial  temperature  distribution  0o(Z*O)  at  K  -  2 
points. 

The  following  numerical  procedure  can  be  used  for  the  solution 
of  the  differential  equation.  For  small  t»  Rs(t)  can  be  expanded 
in  a  Taylor  series  about  x^O* 

•  **  _2  •••  3 

Rs(x)  =  Rs(o)  +  rs(o)t+  r8(o)  jj  +  R#(o)  J7  +  .... 

K 

where  R#(0)  is  obtained  by  satisfying  Eq.  (11). 

.  K-l 

This  can  then  be  used  to  determine  R#(tj)»  1^(xj)»  » •  •  *  Rs(xj) 
while  Eq.  (11)  can  then  be  used  to  determine  R  (x^ ) •  The  nu¬ 
merical  scheme  then  proceeds  as  followst 

t 

Rs(Ti  +  Ax)=  Rs(Ti)  +  AxRs(Ti) 

'R#(Tj  +  Ax)  =  R$(xi  )  +  Ax  R  s(xi  ) 

K-l*  K-l  K 

rs(ti  +  Ax)  =  r  (xt )  +  Ax  R s(xi ) 
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K 

and  finally  Rs(t^  +  At)  can  again  be  evaluate^  from  Eq .  (11).  This 
procedure  is  continued  until  the  entire  solution  is  determined. 

Method  II: 

A  thorough  investigation  of  the  heat  balance  integral  technique 
for  slabs  is  covered  by  Goodman  [2].  Basically,  the  technique  is 
much  the  same  as  the  momentum  integral  of  fluid  dynamics.  The 
heat  conduction  equation  is  satisfied  on  the  average  by  integra¬ 
ting  it  over  the  thickness  of  the  body.  A  second-degree  poly¬ 
nomial  temperature  profile  is  then  assumed  and  the  three  arbitrary 
constants  are  evaluated  from  the  boundary  condition*.  Substitu¬ 
tion  of  tR  is  assumed  profile  into  the  integrated  heat  conduction 
equation  l'eads  for  the  case  of  the  sphere  to  a  second  order, 
non-linear,  ordinary  differential  equation  in  terms  of  the  melt 

radius  r  (t),  and  its  first  and  second  derivatives.  The  solution 
s  • 

is  easily  obtained  by  numerical  procedures. 

Before  integrating  the  heat  conduction  equation  (3),  the  follow¬ 
ing  transformation  is  made: 

Let 


u  =  Tr 


Equation  (3)  then  becomes  (assuming  constant  thermal  properties) 


®u 

9^u 

«t  = 

k  — r 

8r2 

j  b  <  r  <  rg 

• 

k 

where  * 

_  — 

pc 

and  the 

conditions  (4a 

to  e ) 

become 

(12) 
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a)  u(r,tm)  =  uo(r,tm) 


b)  u(rs,t)  =  rsTm 

k  u(rs,t)  dr 

o  0(t)  =  f[(£j)  — r— ]-pi--3r 

S  rs  » t  s 

<>  <Sr> 

Dj  t 

e)  rs(tm)  =  a 

drs 

while  the  condition  (  . .  ) ,  =  0  remains  unchanged.  The  inte~ 

at  t 

gration  of  Eq.  (12)  over  the  thickness  of  the  body  yields 


(13) 


r  n2 
1  s  u 


<-r  s  9u  f  »  i. _ «  , 

J  aTdr  =  V  «  2 

b  b  nr 


Upon  letting: 


r  r  s 

G  =  |  5  u  dr  =  ■  $  Tr  dr 
b  b 


it  can  readily  be  shown  that 


dG  ^_s  r.9u  ,8u.  , 

dt  Tmrs  dt  “  K  A  r  _  .  ^0r^  }  ' 

r  c  t  *  D  »  t 


(14) 


The  substitution  of  Eqs.  (13b»c  and  d)  into  Eq.  (14)  yields 


dfi-T  r 
dt  *m  s  dt 


drs  Kr  dr 

=  —  [Q+ pL— ]  +  k  Tm  -  KT(b,t)  .  (15) 


It  is  now  assumed  that  the  temperature  T(r,t)  can  be  expressed 
at  any  time  as  a  second  order  polynomial  in  (r-b)  with  time 
dependent  coefficients.  That-is» 
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T  =  A  +  B(r  -  b)  +  D(r  -  b)^  . 

Us’nq  conditions  (4b, c  and  d),  A,  B,  and  D  can  be  determined. 
The  resulting  temperature  profile  is 


drs 

Q  +  pL-rr 

T=Tm-C - Sr^Kfr, 


b)  _  Li.  -b)2, 

b)  (r  -b)J 


(16) 


for  rg  4  b 


Therefore,  if  the  solution  for  r  can  be  found,  the  temperature 
profile  through  the  thickness  of  the  body  would  be  specified  at 
any  time  by  Eq.  (16).  G  and  T(b,t)  can  now  be  evaluated  from 
Eq.  (16)  and  the  subsequent  substitution  into  Eq.  (15)  yields 


d2r. 


dt 


dr. 


( r  -  b ) ( 9r  +  7b) 


dr 


Q  Uis,,Vis  -  -r  ,u,  d1t 


( rs  ‘b)  (3rs 


+  12*  (3r  -b) 

- § - }  (17) 

►  5b) 


a  second  order,  non-linear,  ordinary  differential  equation  for 
rs(t). 


The  numerical  procedure  used  here  for  the  solution  of  Eq 
follows  that,  of  Method  I. 


(17) 


For  small  ( t  -  t  ),  r  (t)  is  expanded 
ms 

in  a  Taylor  series  about  t=  t  .  The  first  two  coefficients  of 

IU 

this  series  are  determined  from  the  conditions  at  t=t_,  and  the 

d2rs  m 

third  coefficient  ( ■  r pr ) .  ,  is  obtained  by  evaluating  Eq.  (17) 

dt  d3r  d™r 

The  remaining  coefficients,  (  TTf") *  •••  ("  tf^t-t 


at  t.  =  t 


m* 


dt'' "^t=tm*  dt 

are  obtained  by  successively  differentiating  Eq.  (17)  with 

respect  to  t.  and  evaluating  the  results  at  t  =  t  . 

m 

dr. 

To  determine  rg(t^+At)  and  ( t  ^  +  At )  the  following  numerical 
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scheme  is 


used  t 

rg  ( ti +  At ) 
drs 

—  (ti  +  At) 


=  r s  ( t i )  +  At 
drs 

=  TT(ti)  +  At 


dt‘ 


(ti) 

(ti) 


d2rs 

- r  ( t ,  +  At )  is  determined  from  Eq.  (17). 

dt2  1 


and 
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Numerical  Calculations 


Calculations  using  both  Methods  I  (Citron)  and  II  (Goodman)  are 

performed  for  a  one  inch  thick  7-3/4  inch  outside  diameter  alu- 

* 

minum  sphere  under  a  constant  radial  heat  flux  and  average 

constant  thermal  and  physical  properties  (Tabxa  I).  Two  values 

of  o  (=ha/k=0. 70,  1.10),  the  non-dimensional  surface  heat  trans- 
0 

fer  coefficient,  and  a  stagnation  temperature  of  1850°R,  are  used 
in  the  analysis.  These  values  were  chosen  on  the  basis  of  typi¬ 
cal  heat  inputs  to  re-entry  vehicles. 

The  pre-melt  solution  for  the  temperature  is  obtained  through  the 
use  of  the  one-dimensional  solutions  in  Ref.  [17],  The  outer 
surface  temperature-time  histories  as  well  as  the  times  at  which 
melting  first  occurs  are  shown  in  Fig.  1.  The  corresponding 
temperature  distributions  at  the  melt  time  are  given  in  Fig.  2. 
The  melt  analysis  of  Methods  I  and  II  are  then  used  to  proceed. 

Method  It 

A  six  term  Taylor  series  expansion  for  0(Z,t)  is  assumed.  The 
substitution  of  this  expansion  into  Eq.  (11)  yields  a  non¬ 
linear,  ordinary,  differential  equation  containing  up  to  the 

third  order  time  derivatives  of  R  .  It  should  be  noted  that  the 

•••  s 

highest  order  derivative  (Rg)  appears  linearly  in  this  equation. 
As  previously  mentioned,  this  necessitates  the  matching  of 
0(Z»O)  with  the  pre-melt  solution  (at  t=0)  at  K-2  points.  In 
this  case,  therefore,  only  one  match  point  is  necessary.  The 
point  chosen  is  at  the  insulated  surface  (r=b).  This  enables 

•e« 

the  evaluation  of  Rg(0),  while  Rg(0)  is  found  from  Eq.  (11).  The 
# 

Since  the  surface  temperature  is  assumed  to  remain  at  melt  tem¬ 
perature  during  the  melt  analysis,  a  constant  heat  flux  corres¬ 
ponds  to  steady-state  aerodynamic  heating  during  melting. 
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numerical  procedure  previously  discussed  is  then  used  to  obtain 

Rg  as  a  function  of  x«  Figures  3  and  4  illustrate  the  values  of 

melt  depth  versus  time.  It  should  be  noted  that  the  solution 

dr  s 

for  rg(t)  rapidly  approaches  a  steady  state;  that  is,  ap¬ 

proaches  a  constant  value.  This  is  not  unexpected  since  at 
t=t.m,  the  pre-melt  analysis  yields  a  temperature  profile  through 
the  thickness  which  is  almost  at  the  constant  melt  temperature. 
Therefore,  from  the  boundary  condition  on  the  heat  input  [Eq.  (4c)], 

drs 

Q  — ♦  -  pL  dt  » 


and  since  in  this  analysis  it  is  assumed  Q=  constant. 


Hi 

dt 


— r  ;  a  constant  . 
pL 


This  condition  is  indeed  rapidly  approached  as  shown  in  Figs.  3 

and  4. 


The  Taylor  series  expansion  Eq.  (10)  can  now  be  used  to  determine 
0(Z»x)  versus  Z  for  a  specific  x.  Figs.  5  and  6  show  temperature 
profiles  for  times  greater  than  the  melt  time,  while  Fig.  2  pre¬ 
sents  the  results  at  the  melt  time. 


Method  I I » 

The  numerical  procedure  previously  discussed  is  applied  to  Eq. 
(17).  The  resulting  numerical  values  for  the  melt  depth  are  pre¬ 
sented  in  Figs.  3  and  4.  Equation  (16)  is  then  used  to  calculate 
the  temperature  profiles  at  the  melt  time  (Fig.  2)  and  at  times 
greater  than  the  melt  time  (Figs.  6  and  6). 
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Conclusions 


The  purpose  of  this  report  was  to  obtain  the  one-dimensional 
ablation  characteristics  of  a  hollow  sphere  by  adapting  the 
Taylor  series  expansion  technique  (Method  I)  originally  devel¬ 
oped  for  a  slab  by  Citron  [l]»  and  also  by  using  Goodman’s  [2] 
heat  balance  technique  (Method  II).  Because  of  the  relative 
simplicity  of  the  numerical  procedures  of  Method  II,  it  was 
deemed  appropriate  to  obtain  solutions  using  this  heat  balance 
technique  and  to  compare  these  to  the  more  exact  solutions  de¬ 
termined  from  Method  I. 

A  six-term  Taylor  series  expansion  was  used  in  Method  I.  It 
was  necessary  to  include  this  number  of  terms  in  order  to  ap¬ 
proximate  the  temperature  profile  at  the  melt  time  with  suf¬ 
ficient  accuracy.  A  quadratic  temperature  profile  for  all  time 
was  used  in  Method  II  and  the  comparison  indicates  that  the 
predicted  profile  is  in  good  agreement  with  both  the  exact  pre¬ 
melt  solution  and  the  solution  obtained  from  Method  I  (Fig.  2). 

In  obtaining  a  solution  by  using  Method  I  (six-term  expansion), 
however,  it  was  necessary  to  use  the  insulated  surface  tempera¬ 
ture  at  the  melt  time.  Hence,  the  exact  solution  for  the  tem¬ 
perature  profile  is  approximated  more  closely  by  Method  I  than 
by  Method  II  (Fig.  2).  This  also  accounts  for  the  differences 
between  the  temperature  profiles  shown  in  Figs.  5  and  6. 

The  melt  depths  predicted  by  both  methods,  however,  are  in  ex¬ 
cellent  agreement  (Figs.  3  and  4).  Furthermore,  as  shown  in 
Figs.  3  and  4,  the  rate  of  ablation  rapidly  approaches  a  con¬ 
stant  value.  This  is  not  surprising  since,  as  shown  in  Figs. 

2,  5  and  6,  a  constant  temperature  distribution  is  also  rapidly 
approached. 
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In  conclusion,  it  is  clearly  seen  that  the  numerically  simpler 
technique  (Method  II)  can  be  applied  to  the  ablating  sphere  and 
yields  results  which  are  quite  close  to  those  obtained  using  the 
more  complex  technique  (Method  I).  Both  of  these  methods  pre¬ 
dict  ablation  profiles  which  are  almost  identical.  It  is  appar¬ 
ent  that  future  investigations  should  include  both  the  two- 
dimensional  effects  and  the  effects  of  variation  of  the  material 
thermal  properties  with  temperature. 
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